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Abstract. We investigate theoretically the implementation of two-qubit gates in a 
system of two coupled superconducting qubits. In particular, we analyze two-qubit 
gate operations under the condition that the coupling strength is comparable to or 
even larger than the anharmonicity of the qubits. By numerically solving the time- 
dependent Schrodingcr equation under the assumption of negligible decoherence, we 
obtain the dependence of the two-qubit gate fidelity on the system parameters in the 
case of direct and indirect qubit-qubit coupling. Our numerical results can be used to 
identify the "safe" parameter regime for experimentally implementing two-qubit gates 
with high fidelity in these systems. 
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Two-qubit gate operations in superconducting circuits with strong coupling and weak anharmonicittft. 
1. Introduction 

Superconducting (SC) circuits based on Josephson junctions are promising candidates 
for the realization of scalable quantum computing on a solid-sate platform, due to 
their design flexibility, large-scale integration and controllability (see the reviews in 
Refs. [H El El IH El El [7j). SC qubits, include the charge jH], flux [9], and phase qubits 
[TOj [11] as well as their variants, capacitively shunted flux qubits [12] and capacitively 
shunted charge qubits (transmon) [13]. The phase qubit, the capacitively shunted flux 
qubit and the transmon qubit are relatively insensitive to charge noise and can be 
operated over a wide range of parameters. Single-qubit gates [Ti] . two-qubits gates 
[T5| [To] and simple quantum algorithms [T7] with these types of qubits have been 
demonstrated experimentally in recent years. However, comparing with the flux qubits, 
the common disadvantage of these types of qubits is their weakly-anharmonic energy 
level structure, i.e., the detuning between adjacent transition frequencies is very small. 

Generally, the influence of the small anharmonicity (denoted by A) on quantum 
gate operations can be neglected when the qubit-field or qubit-qubit coupling strength 
is very small compared with A. However, for the practical application of quantum 
computation, one wants to maximize the number of quantum gate operations with a 
given coherence time. In other words, we must implement quantum operations as fast 
as possible, which requires a strong qubit-qubit or qubit-field coupling to be employed 
during the single- and two-qubit gate operations [18J. The anharmonicity of SC qubits 
will influence the quality of quantum gates more and more with increasing coupling 
strength. Recently, there have been a number of theoretical studies analyzing the effects 
of weak anharmonicity of SC qubits on the operation of single-qubit gates and several 
optimization strategies have been proposed based on varying driving pulse shapes and 
sequences [HI [201 Ell E21 [23] . Similar to single-qubit gates, the weak anharmonicity of 
SC qubits will also influence the implementation of two-qubit gates. Then two questions 
arise naturally: (1) how much the weak anharmonicity of the qubits influence the 
implementation of two-qubit gates in a system of coupled SC qubits? (2) how strong can 
the coupling be while allowing a high two-qubit gate fidelity? In other words, how fast 
can two-qubit gates with high fidelity be implemented, given the weak anharmonicity 
of SC qubits? 

Motivated by the above questions, in this paper we study the implementation of 
two-qubit gates with superconducting systems in the strong coupling regime. First, we 
introduce some possible methods for implementing two-qubit gates and qualitatively 
discuss the effect of strong coupling (section II). Then, in section III, we numerically 
simulate the influence of the coupling strength and anharmonicity on the fidelities of two- 
qubit gates in different superconducting systems, and show that the "safe" parameter 
regime for implementing two-qubit gates with high fidelity can be identified, which is 
useful for guiding experimental efforts based on superconducting qubits. Finally, we 
conclude with a brief summary in section IV. 
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Figure 1. (Color online) System with direct (a) and indirect (b) qubit-qubit coupling. 
Here, g, Gj and Aj (j = A, B) are the qubit-qubit, qubit-cavity coupling strength and 
anharmonicity, respectively. 



2. Model and qualitative discussion 

As shown in Fig. 1, as model systems we consider two directly (a) or indirectly (b) 
coupled SC qubits with weakly-anharmonic multilevel structure (such as transmon or 
phase qubits). Here it should be pointed out the flux qubits have a strong anharmonicity, 
and the problem discussed in this paper is not a serious limitation. The two lowest levels 
{\0)j, separated in energy by huj (j = A,B), are the computational basis, and 

the nth (n > 2) higher levels are different from nhujj by he J n . Here e J n has the standard 
nonlinear oscillator form e£ = Aj(n — l)n/2 [21] and Aj is the anharmonicity of the 
qubit, and it is positive in our paper. 

In the case of direct qubit-qubit coupling, two qubits are directly (capacitively) 
coupled, while they are dispersively coupled to a common transmission line resonator in 
the case of indirect qubit-qubit coupling. The Hamiltonian of these two types of coupled 
system is given by (ft = 1) 12H1 123 [2H ESI EOl EH E2! 

N-l 

H direct = J- [(nou A - e£) \n) A (n\ + (nu B - e£) \n) B (n\] + gJ x A ® J&,a) 



n=l 

'N-l 



H indirect = u c a)a + 
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where H d and H ld denote the Hamiltonian for the system with direct and indirect 



Two-qubit gate operations in superconducting circuits with strong coupling and weak anharmonicityA 



qubit-qubit coupling, N is the number of levels in each SC qubit, rf n _ ln = y/n is the 
level- dependent coupling matrix element, and c^_ ln = \n — l)j(n\ + \n — l)j(n\ is the 
effective Pauli spin operators for levels \n — 1) and |n). Also, u c is the frequency of 
the quantized cavity mode; g and Gj denote the qubit-qubit and qubit-cavity coupling 
strength. 

In order to qualitatively analyze the implementation and fidelity of two-qubit gates, 
we assume that each qubit has three levels. Then, the Hamiltonian of direct qubit-qubit 
coupled system ^H direct ) 1 under the rotation- wave approximation (RWA), can be reduced 
to 

Rf rect = J- [ w ,|l),(l| + (2^-A,)|2),(2|] 

j=A,B 

+0[|O1)(1O| + v / 2|02)(ll| + v / 2|20)(ll| + 2|12)(21| + h.c], (2) 

where \mn) denotes \m)A\n)B- 

For the system with indirect qubit-qubit coupling, under the dispersive qubit-cavity- 
coupling condition, i.e., | 5j \ = \ ujj — u c Gj (j = A,B), the qubits will exchange 
energy by virtual photon processes. Then we can obtain the Hamiltonian of the effective 
qubit-qubit interaction by a Frohlich transformation [331 IMl 1351 [36] . 

H indirect = exp ^ S) H id exp(S) 
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Here, we have assumed that Ga = Gb = G. 

The terms proportional to G 2 in the first four terms of equation (3) represent level 
shifts, and the fifth term describes two-photon processes. Under the dispersive qubit- 
cavity-coupling condition, the cavity mode is only virtually excited during the gate 
operation, and therefore the third, fourth, and fifth terms of equation (3) vanish. Then, 
the Hamiltonian (3) can be simplified further as [371 EH ESI EB SI] 

E [^i|l),(l| + (^2-A J )|2) J (2|] 

j=A,B 

[v / 2^ cff ,i|02)(ll| + v / 2p off ,2|20)(ll| + 3cfr,3|01>(10| + 2& ff>4 |12>(21| + /i.e.] .(5) 
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Figure 2. (Color online) The energy-level diagram of two-qubit product states for 
the iSWAP gate (a), and the controlled-Z gate (b) in the system with direct qubit- 
qubit coupling. Red levels denote the states in the computational basis. The black 
dashed arrows are the resonant transitions used for realizing the two-qubit gates and 
the green dotted arrows are the main undesired transitions, which adversely affect the 
implementation of two-qubit gates. The couplings g and \[2g are indicated in blue, 
while the detuning between levels is indicated in black. This figure also applies the 
system with indirect qubit-qubit coupling when the corresponding couplings is replaced 
by ffeff.m (m = 1,2,3,4). 
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Now, we obtain an effective interaction Hamiltonian similar to the Hamiltonian (2) in 
the system with direct qubit-qubit coupling. 

From the Hamiltonians (2) and (5), it is easily seen that various two-qubit gates 
can be realized by appropriately adjusting the qubit frequencies (ua, oob) both in the 
system with direct and indirect qubit-qubit coupling. For example, by setting oja = wb 
{ub = u A + Ab), the resonant transition between state |01) and 1 10) (|11) and 1 02) ) 
can be obtained as shown in Fig. 2. Then the two-qubit iSWAP [T5] (CZ [TBI fT7] ) 
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gate can be realized after an interaction time gt g = n/2 or g e s,zt g = ix/2 {y/2gt = n 
or \/2g c ft t it = 7r). Here it should be pointed out that some undesired transitions [see 
the (green) dotted arrows in Fig. 2] have been neglected in the weak-coupling regime 
g <C \Aj\ or g e s,m *C |Aj| [m = 1 — 4; j = A, B). However with increasing coupling 
strength g or g^m, the average amplitude g/\Aj\ or g eSm /\Aj\ of undesired transitions 
will become larger and larger, which can not be neglected again and will reduce the 
fidelity of the two-qubit gate. So, the relative value of the coupling strength g or g e s t m 
and the anharmonicity Aj is an important parameter for the quality of the two-qubit 
gate. In the two-qubit gate scheme based on SC qubits, a very strong qubit-qubit or 
qubit-cavity coupling strength cannot be employed due to the weak anharmonicity of 
the qubits, if one wants to obtain a high fidelity. How strong the coupling can be, while 
allowing high two-qubit-gate fidelities, will be analyzed in detail in the next section. 

3. Numerical results 

In this section, we will numerically calculate the fidelity of two-qubit gates in the circuits 
with either direct or indirect qubit-qubit coupling. Importantly, the present numerical 
results can help identify the safe parameter regime for implementing two-qubit gates 
with high fidelity. Here, we neglect the noise and decoherence of system in order to 
show explicitly the influence of coupling strength and anharmonicity on the fidelity 
of two-qubit gates. Here, it should also be pointed out that the single-qubit gates 
are performed using microwave pulses (with frequencies of a few of GHz), while the 
frequency tuning for the two-qubit gates are implemented using trapezoidal pulses. 

Here, the fidelity of a two-qubit gate is defined as the Euclidean distance between 
the target Ut and the actual evolution U(t g ) |22j, 

F = l-±\\U T -F*U(t g )P\\l, (7) 

where U(t) is the usual time evolution operator obeying the Schrodinger equation 
U(t) = —j.H(t)U(t) in the full space of the quantum system. Here = tr(X^X) 

where X is an arbitrary operator. P is the projection operator on the two-qubit 
computational {|00), |01), |10), |11)}; 

U T = |00)(00| - z|01)<10| - i|10)<01| + 1 1 1> <1 1 1 

corresponds to the two-qubit iSWAP gate, and 

U T = |00)(00| + |01> (01 1 + |10)(10| - |11)(11| 

corresponds to the two-qubit CZ gate. Here it should be pointed out that single-qubit 
rotations and an overall phase factor = e l0Aa ' z , = e tdBCT ^ , Uj = e ldI are used in 
the numerical calculations in order to eliminate any extra phase factors; I is the unit 
matrix and 

af = |00)(00| + |01)(01| - |10)(10| - |11)(11|, 



Two-qubit gate operations in superconducting circuits with strong coupling and weak anharmonicityl 



of = |00)(00| - |01)(01| + 1 10> <10| - |ii)(ii|. 

Specifically, in our numerical calculations, we replace the unitary operation U (t g ) in 
Eq. (7) by U'(t g ) = UiUfU^U(t g ) and choose 9 A , 9 B and 9 that maximize the fidelity. 

We also note here that in our numerical calculations we do not use the RWA. But, 
there is almost no difference between these results shown below and the numerical results 
with the RWA (not shown in this paper). The reason is that the parameter regime that 
we consider does not reach the ultrastrong coupling regime and thus the RWA is valid 
here. Very recently, the influence of the counter-rotating terms in the Hamiltonian on 
the two-qubit gates in the ultrastrong coupling regime has been studied in a related 
system [12]. Also, the effect of counter-rotating terms were studied in [33] . 

3.1. System with direct qubit-qubit coupling 

In this subsection, based on the original Hamiltonian Eq. (la), we numerically calculate 
the influence of the coupling strength g and anharmonicity Aj on the fidelities of the 
two-qubit iSWAP and CZ gates (see Figs. 3-5). Here we consider the two-qubit iSWAP 
and CZ gates implemented in experiments [15J. In Figs. 3(a) and (b), we plot the 
fidelities of the two-qubit iSWAP gate (FiswAp) and the CZ gate (Fez) as functions of 
g/A B in a circuit with direct qubit-qubit coupling, where we consider each SC qubit 
to have three levels (same approximation will be used in Figs. 4 and 5). From Fig. 3(a) 
and the (green) solid line in Fig. 3(b), it can be seen that the fidelities of these gates 
decrease with increasing g/A B , and the present numerical results can help identify the 
safe parameter regime for realizing two-qubit gates with high fidelities. As shown in 
Fig. 3(a), if we want to implement the two-qubit iSWAP (CZ) gate with fidelity higher 
than 99% (99.2%), the safe parameter regime is g/A B < 0.152 (g/A B < 0.24). In other 
words, based on the relationship gt g = ir/2 for the iSWAP gate and \/2gt g = n for the 
CZ gate, the present numerical results can also identify the time limit for implementing 
two-qubit gates with high fidelity. For example, here the shortest gate time is t g ps 16.4 
ns (t g ~ 12.9 ns) for implementing a two-qubit iSWAP (CZ) gate with fidelity higher 
than 99% (99.2%). 

The (green) solid line in Fig. 3(b) shows small oscillations in the fidelity of the two- 
qubit CZ gate. This result is due to the frequency mismatch between the undesired 
transitions and the resonant transition [see Fig. 2(b)], and it demonstrated that the 
fluctuations of the system parameters will influence the implementation of two-qubit 
gates. Based on the idea of adiabatically eliminating undesired transitions, these 
oscillations can be reduced by slowly adjusting the frequencies of the qubits during 
the gate operation. As shown in the inset of Fig. 3(b), the frequency of qubit B starts at 
l.lco B , is first ramped down to u B in r d , then ramped up to l.lu B after an interaction 
time t g (\/2gt g = n). During the full gate operation time (2t<2+£ s ), the frequency of 
qubit A is fixed. Using such pulses, we numerically calculate the fidelities of the two- 
qubit CZ gate for different values of r d and present the results in Fig. 3(b) [See dashed, 
dotted and dot-dashed lines in Fig. 3(b)]. It can be seen that the oscillations of the 
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Figure 3. (Color online) The fidelities of the two-qubit iSWAP (a) and CZ (b) gate as 
functions of g/ Ab in a circuit with direct qubit-qubit coupling. Some representative 
dots are denoted by the dashed lines and red circles in order to present the relationship 
between the gate time t g and fidelity F. The red arrows point out the parameter regime 
corresponding to two-qubit gate with high fidelity. In figure (b), the qubit frequencies 
are adiabatically adjusted during the gate operation, as shown in the inset part. The 
system parameters used here are: (a) loa/2it = 5.5 GHz, lob — wa, Aa/2-k = 0.15 GHz, 
and A b /2tt = 0.1 GHz; (b) u} A /2n = 7.16 GHz, A A /2n = 0.087 GHz, A b /2tt = 0.114 
GHz, and ujb — + Ab- 
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Figure 4. (Color online) The fidelities of the two-qubit iSWAP gate (a) -FIswap and 
CZ gate Fqz (b) versus A^/g and A.s/g in a circuit with direct qubit-qubit coupling. 
The dashed lines correspond to the parameter regime for implementing a two-qubit 
gate with fidelities 95% and 99%. The system parameters are the same as in Fig. 3 
except for g/2ir = 0.2 GHz. 




Figure 5. (Color online) The fidelities of the two-qubit iSWAP gates versus A^ and 
g in a circuit with direct qubit-qubit coupling. The dashed lines correspond to the 
parameter regime for implementing two-qubit gate with fidelities 95% and 99.5%. The 
system parameters are the same as in Fig. 3 except for A^ = Ag. 



fidelity can be eliminated by adiabatically adjusting the qubit frequencies during the 
gate operation. This numerical result provides a method to reduce the influence of 
parameter fluctuations on the implementation of two-qubit gates. 
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Figure 6. (Color online) The fidelity of the two-qubit CZ gate versus ^a/GsSIj 
Ag /<7 ff i (a) and versus g c g j , (b) in the system with indirect qubit-qubit coupling. 
The dashed lines correspond to the parameter regime for implementing two-qubit gate 
with fidelities 95% and 99%. The basal system parameters are: uj c /2tt = 6.9 GHz, 
lu a /2it = 8.2 GHz, lu b = lu a + A B , 5j = u d - lo c (j = A, B); And G = 0.2 GHz for 
panel (a), Aa/2tt = A B /2n GHz for panel (b). 



In order to show the influence of A^ and A B on the two-qubit gates, we plot the 
fidelities of the two-qubit iSWAP and CZ gates as functions of A a/ g and A B /g in Fig. 4. 
It is easily seen from Fig. 4(a) that the anharmonicities A^ and A B have equal effects 
on the two-qubit iSWAP gate, i.e., the larger the anharmonicities Aj- (j = A,B) are, 
the higher the fidelity. This symmetric property disappears in the two-qubit CZ gate 
due to the asymmetry in the condition on the parameters, ub = uja + A b [see Fig. 4(b)]. 
In other words, the influence of the anharmonicity A^ on the two-qubit CZ gate can be 
neglected when w B = ua + A b is chosen. In addition, the dashed lines in Fig. 4 indicate 
the safe regime of Aj/g (j = A,B) for implementing two-qubit iSWAP and CZ gates 
with fidelity higher than 99%. 

In Figs. 3 and 4, either the anharmonicity Aj or the coupling strength g have been 
set to a fixed value. A natural question is whether the conclusions obtained from Figs. 3 
and 4 are universal. In other words, will the properties of Figs. 3 and 4 change much 
when either Aj or g is changed? Thus, we now present in Fig. 5 three-dimensional (3D) 
plots of the dependence of Inswap on g and A#. It is shown that the fidelity of two-qubit 
gates are approximately determined by the ratio of the qubit-qubit coupling strength g 
to the anharmonicity Aj of the SC qubits. As a result, the conclusion obtained from 
Fig. 3(a) [or Fig. 4(a)] will not be changed when adjusting A B (or g). A similar property 
is also obtained from the two-qubit CZ gate (the corresponding figures are not shown 
in this paper because are very similar to Fig. 5). 
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Figure 7. (Color online) The fidelities of the two-qubit gates as a function of <?/A_g 
(a) and g e s,x/^-B (b) in systems with direct (a) and indirect (b) qubit-qubit coupling, 
when the three, four, or five lowest levels are considered for each qubit. The system 
parameters are the same as in Fig. 3 or 6. The green, red circles in (a) and cyan circle 
in (b) mark respectively the experimental parameters regime in Refs. |15j . |16j . |17j . 



3.2. System with indirect qubit-qubit coupling 

In this subsection, based on the Hamiltonian Eq. (lb), we present the results of numerical 
calculations for the dependence of the fidelity of the two-qubit gates on the effective 
qubit-qubit coupling g effl and anharmonicity Aj of SC qubits. Here the two-qubit CZ 
gates are realized based on the qubit-cavity dispersive interaction method [XT] , and the 
parameter 

91 ( 1 ± \ - 

geStl 2 \5 B - A B 8 A ) 5 a 
under the condition = + A^. 

In Fig. 6, we present the 3D plots of the dependence of F C z on AA/g e s,i and A B /g e ff,i 
[panel (a)], and g> e ff,i and A B [panel (b)], where we consider the SC qubits to have three 
levels. Using dashed lines, we have denoted the parameter regime for implementing 
two-qubit CZ gate with fidelities 95% and 99%. It is shown from Figs. 6(a) and (b) that 
high-fidelity areas correspond to the weak-coupling regime g e s,i/ Aj <C 1 (j = A,B), 
while low fidelity corresponds to the strong-coupling regime, where g e ff,i is comparable to 
or larger than Aj . This property is similar as that in the system with direct qubit-qubit 
coupling. The present numerical results can be used to identify the safe parameter 
regime for implementing the two-qubit CZ gate with high fidelity in the circuit with 
indirect qubit-qubit coupling. 

3.3. Going beyond the three-level approximation 

Until now, three-level-system approximation for qubits has been used in the above 
numerical calculations. It is then natural to ask the following question: will our 
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conclusions, obtained from the above numerical results, still be valid for qubits with 
iV (N >3) levels? To explore this, in Fig. 7, we plot the fidelities of the two-qubit 
iSWAP and CZ gates as functions of g/A B (or g c s^/A B ) in the system with direct (or 
indirect) qubit-qubit coupling when each qubit has three, four or five levels. It can be 
seen from Fig. 7 that there is not much difference between the numerical results based 
on the three-, four- and five-level approximations for the qubits. So, our conclusions 
obtained from the above numerical calculations are still valid for iV-level (with N >3) 
SC qubits. 

3.4- Limits on the gate fidelities of recent experiments imposed by weak anharmonicity 

In order to serve as a guide for future experiments, we compare our numerical results 
with corresponding experiments and show the limited fidelity of two-qubit gate based 
on SC qubits with weak anharmonicity. Based on the experimental parameters (ua/27t, 
co B /2n, A A /2n, A b /2tt, gj2ix) equal to (5.5, 5.5, 0.15, 0.1, 0.011) GHz and (7.16, 7.274, 
0.087, 0.114, 0.0091) GHz, two-qubit iSWAP pj and CZ [16] gates with fidelities 63% 
and 70% were implemented in the circuit with direct qubit-qubit coupling. In the circuit 
with indirect qubit-qubit coupling, a two-qubit gate [17] with fidelity 85% was realized 
with system parameters (u c /2n, ua/2tt, u b /2tt, A^/27r, A B /2n, Ga/2ti = G B /27i) 
equal to (6.9, 8.2, 8.45, 0.2, 0.25, 0.199) GHz. Corresponding to the above experimental 
parameters, in Fig. 7 we indicate the ideal fidelity (see the green, red and magenta circles) 
based on our theoretical calculations. From the comparison between experiments and 
our numerical calculations, we show that two-qubit gates with fidelities 99.52%, 99.91%, 
and 99.2% can be realized, in principle, if the influence of decoherence can be eliminated. 
Recently, the effects of decoherence on quantum gates and possible optimization routes 
were also studied in Ref. |4*4"] . 

4. Conclusion 

We have studied the performance of two-qubit gates in a system of two coupled SC 
qubits under the condition that the coupling strength is comparable to or larger than 
the anharmonicity of the qubits. First of all, by using the three-level approximation 
for the qubits, we analyzed and numerically calculated the dependence of the two- 
qubit gate fidelity on the qubit-qubit coupling strength and the anharmonicity of the 
qubits. Based on extensive numerical results, the safe parameter regime was identified 
for experimentally implementing two-qubit gates with high fidelity. Secondly, we 
numerically calculated the fidelity of the two-qubit gates in the case of four- and five-level 
approximations for the qubits, and demonstrated the validity of our numerical results 
for N- level qubits with N > 3. Our results can serve as a guide for future experiments 
based on SC qubits. 
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